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PREFACE 


The purpose of this book is to review the present knowledge of the 
mixed state of superconductivity. It is the hope of the author that this 
review will be helpful to persons unfamiliar with the subject. The 
fundamentals are discussed in detail to provide sufficient background 
for the understanding of more recent developments. The emphasis is 
on the experimental and theoretical aspects of flux creep and flux 
flow. The thermal properties of the mixed state have not been discussed. 

This publication was originally prepared in partial fulfillment of the 
requirements for a Master of Science degree at John Carroll University, 
Cleveland, Ohio. 
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INTRODUCTION 


The discussion begins with a review of the fundamental properties of 
superconductivity. Then the Ginzburg-Landau equations are described 
and Abrikosov’s solution is discussed. The vortex state is compared to a 
laminar state to determine which state is energetically more favorable. 
The experimental evidence that substantiates the vortex state is given. 
The results of tube magnetization experiments are presented, and the 
flux creep theory is derived. The flux flow models of Bardeen-Stephen 
and Nozieres-Vinen are compared with each other and with the experi- 
mental evidence. Three models for the viscosity coefficients are pre- 
sented. Finally, the rigidity of the Abrikosov structure and the possible 
flux line pinning mechanisms are mentioned. 


Zero Resistance Effect 

In 1908 Kamerlingh Onnes liquefied helium for the first time. Elec- 
trical resistivity was not very well understood and Onnes decided 
to make use of his liquid helium capability to study the dependence, of 
resistivity on temperature. He believed that if he could obtain a metal 
in a pure enough condition its resistivity would approach zero at absolute 
zero. For this reason he chose mercury. Figure 1 shows the results 
obtained by Onnes. 

Onnes found that the resistivity dropped very abruptly to an un- 
measurably small value at a temperature just slightly above the normal 
boiling point of helium (ref. 1). This transition temperature is now 
called the critical temperature TV Initially Onnes explained that this 
result occurred because of the purity of the mercury. He discarded 
this conclusion when he found that the resistivity of mercury disappeared 
even when impurities were added. The addition of impurities did, how- 
ever, increase the temperature span over which the resistivity dropped 
to zero. The sharpness of the transition and its relative insensitivity 
to alloying forced Onnes to conclude that mercury had undergone a very 
fundamental change to a new electronic state. He named the new state 
the superconducting state. 

The discovery of the superconducting state opened a whole new area 
of research. One obvious question was how many other substances 
exhibit the zero resistance effect. As of 1963, more than 900 super- 
conducting materials have been found (ref. 2). Among the known super- 
conducting materials are metals, semiconductors, alloys, and compounds. 
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FIGURE 1. — Resistance of mercury plotted against absolute temperature (ref. 105). 

The new electronic state was also probed experimentally in many 
ways. Onnes (refs. 3 and 4) discovered that if an electric current through 
a superconductor exceeds a critical value the material becomes resistive. 
He also discovered (ref. 5) that if a magnetic field is applied to a super- 
conductor the resistance returns when the field exceeds a particular 
critical value H c . The detailed nature of the magnetic transition became 
known somewhat later (refs. 6 and 7). The critical magnetic field was 
found to be temperature dependent. In many materials the temperature 
dependence is described fairly well by 


Hc = H 0 



( 1 ) 


where H 0 is the critical field at absolute zero. Some experimentally de- 
termined critical field against temperature curves may be found in ap- 
pendix A. In the region above the critical curve the material is in the 
normal state, and in the region below the curve it is in the supercon- 
ducting state. It was suggested by Silsbee (ref. 8) that the critical current 
and the critical magnetic field are related. Silsbee’s hypothesis is that 
the critical current is the current which generates a self field equal to 
the critical magnetic field. 

The superconducting state is very sensitive to impurities and de- 
formations. Figure 2 shows the changes that occur in the resistive 
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Pure single crystal tin 

—Pure polycrystal tin 

Less pure polycrystal tin Current, 



(a) Influence of physical properties of (b) Influence of measuring 

the samples. current. 

Figure 2. — Reduced resistance against absolute temperature (ref. 102). 


transition which depend on the sample condition and the measuring 
current. These measurements, which were made by deHaas and Voogd 
(ref. 9), showed that the transition could take place in a temperature 
span as small as 0.001 K. The specific heat of a superconductor changes 
discontinuously at the critical temperature, and the thermoelectric 
power of superconductors is equal to zero. 

Since superconductivity was discovered, various researchers have 
attempted to measure a resistance in the superconducting state. The 
early potentiometric measurements of Onnes were accurate enough to 
show that resistivity in the superconducting state was 10~ 6 less than 
the normal state resistivity. Since then various investigators have been 
able to put an upper bound on the resistivity by inducing a current in 
a ring of superconducting material and measuring its decay. Onnes was 
the first to use the method. Grassman (ref. 10) reported a resistivity of 
less than 10~ 12 of the normal state value. Recently File and Mills (ref. 11) 
have performed experiments where the decay time of a niobium zir- 
conium (NbZr) solenoid was measured. They found that the decay time 
of the persistent currents must be greater than 100 000 years. 


Meissner Effect 

In 1933, Meissner and Ochsenfeld (ref. 12) discovered an abrupt 
change in the apparent magnetization of tin as it was cooled below the 
transition temperature. The original purpose of their experiment was 
to determine if the current in a superconductor flowed on the surface 
or in the bulk. The experimental scheme was to measure the magnetic 
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field in the vicinity of a U-shaped current-carrying wire. They hoped 
to be able to distinguish between the two modes of current flow by meas- 
uring small differences in the magnetic field. What they found instead 
was that the magnetic field surrounding the wire increased sharply 
when their sample became superconducting. 

Later experiments showed that the magnetic field is totally expelled 
from the bulk of a superconductor when the magnetic field is less than 
the critical field. The applied magnetic field H a is excluded by super- 
conducting surface currents. If these currents are represented by an 
equivalent magnetization, the Meissner effect can be shown graphically 
by plotting magnetization as a function of applied magnetic field (fig. 3). 
Below H c there is no magnetic field inside the superconductor and the 
magnetic flux density B equals zero. Therefore J? = /f a + 47rM becomes 
Ha~ — ^TrM. Above H c the superconductor can no longer expel the 
magnetic field and the sample undergoes a transition to the normal 
state. 



Applied magnetic field, H a 

Figure 3. — Ideal magnetization against applied magnetic field for superconductor. Below 
H c superconductor is perfectly diamagnetic, H„=—4arM. 


The Meissner effect may be complicated by the sample geometry 
and the orientation of the magnetic field with respect to the sample. 
These effects, called demagnetization effects, are discussed in detail 
in references 13 and 14. Basically, figure 3 would be obtained only for an 
infinitely long cylindrical sample with the magnetic field paralled to its 
axis. In this geometry the magnetic field is paralled to the internal mag- 
netization. The critical field H c is defined for this geometry. Any other 
geometry exhibits a lower value of critical field. When the sample 
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geometry or the sample orientation with respect to the magnetic field 
is such that magnetic field penetrates the sample below H c , the sample 
is said to be in the intermediate state. The intermediate state has been 
optically observed using an experimental technique that takes advantage 
of the Faraday rotation effect (ref. 15). Some pictures taken by Faber 
(ref. 16) are shown in appendix B. These particular pictures are of an 
aluminium plate. The dark areas are superconducting regions. The 
pictures show that the intermediate state is composed of alternating 
normal and superconducting regions which are laminar-like in appear- 
ance. The width of the superconducting laminae are seen to decrease 
with magnetic field. The size of the normal and superconducting laminae 
varies because the total magnetic flux through the sample must be 
equal to the applied flux, and the magnetic field in the normal regions is 
H c . 


London Equation and Penetration Depth 

Prior to the discovery of the Meissner effect, Becker, Heller, and Sauter 
(ref. 17) had analyzed the electric and magnetic properties of a “perfect 
conductor” (zero resistance, but no Meissner effect). Their analysis 
predicts that any magnetic field in a material when it becomes super- 
conducting will be trapped inside the sample. This prediction was 
contradicted by the Meissner effect. 

F. London and H. London (ref. 18) postulated a fundamental con- 
tradiction between the Becker-Heller-Sauter prediction and the Meissner 
effect. They proposed that in superconductors 

(2) 

me 

where J is the current density, n s the number of superconducting 
electrons, e the charge of an electron, m the mass of an electron, c the 
speed of light, and A the magnetic vector potential. Then, starting with 
the Maxwell equation 

Vxfl=^/ (3) 

£> 

and taking the curl of both sides of both the London and the Maxwell 
equations yields 

Vxj=-^(Vxi)=-^H 

me me 


(4) 
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and 

V x VXff=~(V x/) (5) 

Substituting equation (4) into equation (5) results in 

vx Vx/r=— 4 tt^H (6) 

me 1 

For a magnetic field parallel to the surface of a superconducting half 
space, equation (6) is reduced to one dimension and has solutions of the 
form 

H=H a e-* lkL (7) 

where k L = (mc 2 /47rn s e 2 ) 1/2 and H a is the applied magnetic field. The 
x coordinate is perpendicular to the surface of the superconducting 
region. This result agrees with the Meissner effect because it predicts 
that in the bulk the magnetic field is essentially zero. The parameter 
X/„ which is called the London penetration depth, is the distance from 
the surface of a semi-infinite slab to the point where H = HJe . The 
shielding currents mentioned earlier reside primarily within a few 
penetration depths of a sample’s surface. A typical value for Xl is 
about 500 A at 0 K. The penetration depth is relatively constant except 
near the critical temperature where n s goes to zero. Table I (p. 65) con- 
tains measurements of the penetration depth for several materials. 


Mixed State 

In 1934 and 1935 Mendelssohn and his coworkers at Oxford observed 
some rather disconcerting experimental results on superconducting 
alloys. They were investigating the Meissner effect in various materials: 
mercury (Hg), lead (Pb), tin (Sn), tantalum (Ta), tin bismuth (SnBi), lead 
bismuth (PbBi), and lead telluride (PbTL). They sought to observe the 
sharp and reversible changes in the specific heat of superconductors, 
which were predicted by the thermodynamic treatments of Gorter and 
Casmir (refs. 19 and 20). Gorter and Casmir assumed that the Gibbs 
free energy of a superconductor G s was related to the free energy of a 
normal conductor G n by 


M 2 

Gn = G*( 0)+g 


( 8 ) 


and had a magnetic field dependence of 
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G S {H) = G S ( 0)+— (9) 

From these assumptions they derived a discontinuity in the specific 
heat. Keeley and Mendelssohn (ref. 21) found that only mercury had a 
reversible magnetization. Tin and lead exhibited some hysteresis; 
tantalum was found to have even more. In addition, the alloys were de- 
scribed by Mendelssohn as being completely irreversible. There were 
other differences between the experimental results for alloys and the 
same results on elements. The empirical fact observed for elements, that 
the resistance and magnetic induction go to zero at the same magnetic 
field value, was not found to be true for alloys, Gorter’s prediction of a 
specific heat discontinuity had been observed by Mendelssohn and Moore 
in tin (ref. 22) and by Keesom and Kok (ref. 23) in thallium. Their specific 
heat measurements were consistent with each other and with Gorter’s 
prediction. But when Mendelssohn and Moore (ref. 24) investigated the 
specific heat of PbTh, which they estimated should have a large specific 
heat discontinuity, they found the result shown in figure 4. No indication 
of a discontinuity was observed. These results strongly suggested that 



FIGURE 4. — Specific heat C p of PbTl 2 against absolute temperature (ref. 24). 
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a fundamental difference exists between the superconducting properties 
of pure metals and alloys. 

It has since been shown (ref. 25 ) that even if the magnetization as a 
function of applied magnetic field curves of alloys are completely re- 
versible they are fundamentally different from figure 3 . Figure 5 shows 
the magnetization that characterizes most alloys. Up to H c i the magnetic 
field is completely excluded, but between H c \ and H c 2 the magnetic 
field inside the material increases until the material becomes completely 
normal at H c 2. Materials which exhibit this type of magnetization curve 
are called type II superconductors. Superconductors that have a mag- 
netization curve as shown in figure 3 are called type I. When the mag- 
netic field of a type II superconductor is between H c 1 and H c 2 the super- 
conductor is said to be in the mixed state. The mixed state is distin- 
guished from the intermediate state in that it is not a geometry effect, 
the mixed state is rather an independent equilibrium state of the 
system. 



Applied magnetic field, H a 

Figure 5. — Magnetization against applied magnetic field for type II superconductor. 


Figure 6 shows how the magnetization is affected by alloying. Starting 
with pure lead, which is a type I superconductor, increasing amounts 
of indium are added. The resulting alloy is a type II superconductor. 
As the percent of indium in the alloy is increased, H c 1 decreases and H c 2 
increases. The table accompanying the figure shows how the ratio of 
H C 2 to Hd varies. 
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Applied magnetic field, H a , tesla 


Figure 6. — Effect of alloying on magnetization curve of lead (ref, 105). Key shows variation 
of Hcu #c2, k , and HcilHci with alloying; thermodynamic critical field, = 0.0550 
tesla. 


Bardeen-Cooper-Schrieffer Theory and Coherence Length 

The modern era of superconductivity began in 1950 when Frohlich 
(ref. 26) suggested that superconductivity might be explained by phonon 
coupling between electrons. Frohlich was able to predict that the de- 
pendence of the critical temperature on the isotropic mass was T c 01 
About the same time, this relation was discovered experimentally in 
mercury (refs. 27 and 28). These results strongly suggested that super- 
conductivity results from an electron-lattice interaction. In 1957, Bardeen, 
Cooper, and Schrieffer (B.C.S.) put forth a theory based on the electron- 
phonon interaction suggested by Frohlich (refs. 29 and 30). The B.C.S. 
theory proposes that the interaction occurs between electrons with oppo- 
site momenta and spin. 

The interaction may be visualized as follows. An electron moves 
through the lattice attracting the positive lattice toward it. A second 
electron can be attracted to the deformed region by the excess of posi- 
tive charge there. This effectively attractive interaction between the 
two electrons results in an energy gap between an electron in the B.C.S. 
ground state and the Fermi level. The energy gap A, which is related to 
the critical temperature by A 01 ksTcy leads to the critical fields and the 
thermal properties mentioned earlier. The attractive force between 
the electrons extends over a distance which is called the coherence 
length The coherence length can also be interpreted as the distance 
over which the energy gap A cannot change drastically. The B.C.S. 
theory results in the following fundamental relation at T= 0 K between 
the energy gap and the coherence length £ 0 : 
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where V/ is the Fermi velocity and h is Planck’s constant divided by 
27 r. Some typical values of the coherence length are exhibited in ap- 
pendix C. Prior to the formulation of the B.C.S. theory, Pippard (ref. 31) 
had predicted that a characteristic length like the coherence length 
must exist. He further proposed that the coherence length is related to 
the mean free path A by the equation 


i=i+i 

£ £. a 


( 11 ) 


Two limitations of the B.C.S. theory are first that it does not predict 
the mixed state of type II superconductors and second that it does not 
apply to current carrying superconductors. 


Abrikosov's Analysis of the Mixed State 

Abrikosov (ref. 32) showed that in the mixed state the magnetic field 
penetrates the superconductor in discrete flux tubes which are enclosed 
in vortices of superconducting electrons. Each tube of flux contains 
a flux quantum <p 0 equal to 2.07 X 10" 15 meter squared tesla. The exist- 
ance of the flux quantum was originally predicted by London (ref. 33, 
p. 151; see also ref. 13, p. 32). Figure 7 shows the spacial variations 
of the current density of the vortex |J|, the magnetic field H , and the 
energy gap A in the region near a tube. Abrikosov’s analysis gives the 
following relations between H c i, H c 2 , and H c : 


and 



(12) 

(13) 


These results suggest that the relationship between the penetration 
depth X and the coherence length £ determines whether a superconduc- 
tor is type I or II. If X/f > 1/V2 then equation (13) implies that H C 2 >H C 
and the superconductor must be type II. Equation (12) predicts that for 
a type II superconductor, H cl <H c . In figure 6, H cl decreased with 
alloying and H c 2 increased. This figure can now be understood in terms 
of equations (11) to (13). Alloying lead with indium reduces the mean 
free path and equation (11) predicts that £ will also be reduced. If f is 
reduced while X remains constant, then equations (12) and (13) predict 
the experimental observation in figure 6. The ratio of X to £ is called the 
Ginzburg-Landau kappa k. 
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r 


r 


Figure 7. — Spatial variations of \J | , H , and A that make up an Abrikosov vortex. Energy 

gap A is proportional to n s . 


GINZBURG-LANDAU THEORY 

The Ginzburg-Landau (G.L.) theory (ref. 34) is a thermodynamic 
approach to superconductivity. Ginzburg and Landau postulated that an 
additional term needed to be added to the Gibbs free energy G S (H) 
(given by eq. (9)). The new term is needed to incorporate the kinetic 
energy (V^r) 2 that results when n s vary spatially. In the Ginzburg- 
Landau formulation, iff is called the order parameter, but it is now 
known to correspond to the wave function of a pair of superconducting 
electrons. The number density of superconducting electrons n s is pro- 
portioned to To keep their theory gage-invariant, Ginzburg and 
Landau included a gradient term and a constant term. The total Gibbs 
free energy is 


c,(«)-c.(o )+ f + -L 




(14) 


The G s ( 0) term can be written as a function of if) utilizing the Landau- 
Lifshitz general theory of second-order phase transitions (ref. 35). 
The Landau-Lifshitz theory states that near the transition temperature 
the Gibbs free energy can be expressed as a series expansion in |^| 2 . 
Therefore, 

G s (0) = G n (0) + a(T)M 2 +^W + - . . (15) 


The minimum value of G s (0) corresponds to the equilibrium value of 
\.iff\ 2 . Both a(T) and fi(T) are temperature-dependent coefficients. 
The minimum of G s ( 0) is obtained by taking the partial derivative of 
£ s (0) with respect to \iff\ 2 and setting it equal to zero: 




w=- 


*(T) 

W) 


(16) 

(17) 
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Equation (18) is obtained by substituting |i/i| 2 into equation (15) and 
comparing the result with equation (8): 


a-(T) m 

2)8(7’) 8tt 

Near T c , if it is assumed that 
and 

P(T)=P(T c ) =/3c 


(18) 


(19) 


then substituting for a{T) and jS(T’) in equation (18) results in 

tt , _ 47T (7’c -70 2 (da\ 2 /o m 

Hc ~Tc ( a7 ) T = Tc (20) 

This is consistent with experimental results. Near Tc, H c is usually found 
to be proportional to T c — T. This agreement tends to validate the as- 
sumptions in equation (19), 

The Ginzburg-Landau equations can be derived by integrating equation 
(14) over the volume of the superconductor and then minimizing the 
result with respect to both xjj and A . This can be done using the Euler 
method (ref. 36). The general results (ref. 34) are 


and 


(?-) 




= \M 2 A + - 4/Vt If*) 


( 21 ) 

(22) 


Both equations (21) and (22) have been reduced to a standard dimension- 
less form where lengths are normalized to X, magnetic field is normalized 
to V2 H c , and i /j is normalized to \\f H = o- The Ginzburg-Landau parameter 
k equals V2 (elhc)H c k 2 . London’s equation (eq. (2)) is obtained from 
equation (22) in the limit of low magnetic field, as H — > 0, because 
i and Vt^=0. 

The Ginzburg-Landau equations have been derived from the micro- 
scopic theory of Bardeen, Cooper, and Schrieffer by Gorkov (ref. 37). 
These equations can be classified as coupled nonlinear equations in iJj 
and A . Because of their inherent complexity they have been solved 
primarily for one-dimensional problems. A partial listing of the problems 
to which they have been applied includes the properties of super- 
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conductors at high fields, the critical fields of thin films, supercooling 
in type I superconductors, the dependence of penetration depth on field, 
and the calculation of the interphase boundary energy. 

It is informative at this point to derive the Ginzburg-Landau equations 
for a simple one-dimensional model. First, insight is gained into the 
Ginzburg-Landau equations. Second, the very important prediction 
that there are two distinct kinds of superconductors can be readily 
derived from the solutions. The pedagogical development which follows 
is in large part due to Newhouse (ref. 14). Consider an infinite super- 
conducting half space x < 0 (fig. 8). If H is in the ^-direction and ,4 is in the 
y-direction, then A ~A(x) and i// = i//(a;). Therefore, H—dA/dx . 
Because Vtff is in the ^-direction, /4*Vi// = 0. Incorporating these facts 
into the total Gibbs free energy equation (eq. (14)) yields 

cm - CM + 1 (g )■ + £ [*<f )’ + (^ +)’] (23) 

The corresponding Euler equations (ref. 36) are 


BGs(H) d / dG s (H) \ 
dijj dx \ dt/Tr / 

and 

dG s (H) d ( dG,(H) \ 
SA dx \ BA X ) 


(24) 


(25) 


where tyx-diff/dx and A x = dA/dx. When 4s A, and A x are regarded 
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as independent variables, equations (24) and (25) can be evaluated 
by taking the appropriate derivatives of equation (23). The result is the 
Ginzburg-Landau equations for this geometry: 


dty m dG s ( 0) 

dx 2 h 2 dijj + \hc)' 1 ' 


(26) 


and 


d-A I6rre 2 

-y~r= — ~ 
ax- mo- 


rn 


To completely determine equation (26), d£ s (0)/di// must be evalu- 
ated exactly. The oc{T) and fi(T) in equation (15) can be expressed in 
terms of experimentally observable parameters. Starting with the 
Maxwell equation. 


V'H— — J 

c 


(28) 


If V • A = 0 then V X VX A— — V 2 A = — drA/dx 2 in one dimension, and 


d 2 A __ 4 77 

dx 2 c 


(29) 


Substituting the result into equation (27) yields 


J = - 


4ef 

me 


\fj 2 A 


(30) 


Equation (31) is obtained by substituting equation (17) into equation (30) 
and comparing with the London equation (eq. (2)): 


oc(T) m_ / c \ 2 

P(Tf~ 4ir \ 2eXL/ 


.(31) 


The penetration depth k L — (mc 2 /47m s e 2 ) 1/2 was defined by equation (7). 
For convenience in writing and reading the equations that follow, 
T= (2 ek L /c) 2 . Expressions for a(T) and fi(T) are obtained by combining 
equation (31) with equation (18) and solving for a{T) and (3(T): 
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and 


a(r>- 


-Sr 

m 


J8(T) 


47 tH 2 
m 2 


r 2 


(32) 


Substituting into the Gibbs free energy equation (eq. (15)) yields 

( 1/2 \ / o _IJ2 \ 

___£ r )^+(-^p)^ (33) 

Taking the derivative gives 


dG g (0) 

dip 


dG n 

"ai/T 


2//| 

m 


15//-}- 


8rrH 2 

Til 1 


r 


(34) 


If dG w (0)/cJx— 0 is substituted into equation (34) and the result is sub- 
stituted into equation (26), 


where 


and 



A 2 



2V2e// c X| 

K ~ he 


(35) 


(36) 


Equation (35) can be put in the same form as the Schroedinger equation 
for a simple harmonic oscillator with real periodic solutions if H ^ // c 
and ip 0 > ip: 


d 2 \p __ k 2 A H 2 x 2 \ 

A! V 


(37) 


The solutions are real and periodic only if 


H = kV2 

H c 2 n + 1 


n = 0, 1, 2, . . . 


(38) 


The maximum value of HjH c occurs when n = 0. Therefore, the lowest 
magnetic field a superconducting material can be in without the nuclea- 
tion of superconductivity is 
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H c2 = kV2H c (39) 

From this equation it can be seen that at k > 1/V2 the upper critical 
field H c 2 becomes greater than H c . This range of k corresponds to type 
II superconductivity. The value k = 1/V2 is then the dividing point be- 
tween types I and II superconductors. 


MIXED STATE, VORTEX OR LAMINAR? 

The theoretical confirmation of type II superconductors from the 
Ginzburg-Landau theory raises the question of the nature of the mixed 
state. What form does the flux inside the superconductor take? Histori- 
cally two models were proposed, the laminar model (ref. 25) and the 
vortex model (ref. 32). In both models the mixed state is thought of as 
consisting of normal regions, where most of the flux resides, and super- 
conducting regions, where London’s equation is obeyed. In the laminar 
model the normal regions are plane laminae of thickness 2£. The vortex 
model proposes that the normal regions, often called normal cores, are 
cylinders of radius The laminar model was proposed because laminae 
have been observed experimentally (appendix B) in the intermediate 
state of type I superconductors. The vortex model resulted from 
Abrikosov’s solution of the Ginzburg-Landau equations. 


Abrikosov Solution of Ginzburg-Landau Equations 

The most important solution of the Ginzburg-Landau equations from 
the standpoint of providing insight into the mixed state was found by 
Abrikosov (ref. 32). Onsager (ref. 38) and Feynman (ref. 39) predicted 
that vortices should exist in superfluid helium. Abrikosov became con- 
vinced that vortices of current must exist in a superconductor. He 
therefore set out to find a doubly periodic solution to the Ginzburg-Landau 
equations. He started out with very restrictive assumptions. He assumed 
that the superconductor extends throughout all space. The supercon- 
ductor is assumed to be near T c , which implies that \ip\ 2 1. The mag- 

netic field inside the superconductor H is set equal to the applied field 
H a . It is assumed that H n is just slightly less than H c 2 . It is further 
assumed that k > 1/V2. 

Abrikosov first sought to solve equation (21). If H a is in the z-direction 
and A is directed along the y-axis, then 


H = Ha 
A = HaX 


(40) 
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Substituting this result for A into equation (21), ignoring the term, 
and taking i jj — ipix) give 


dx 2 


k 2 (-H 2 x 2 + 1)xIi = 0 


( 41 ) 


This equation has solutions of the form 

ifj = e 2 


( 42 ) 


This simple analysis lead Abrikosov to a less restrictive solution of 
equation (21): 


¥ = S C n e ikn ^nM 

71 — — 00 


( 43 ) 


where k and C„ are arbitrary constants. Substituting this solution into 
equation (22) and solving for A yield 


A=H,,x~y k J M 2 dx 


( 44 ) 



FIGURE 9. — Plot shows that for Abrikosov’s solution has the symmetry of a square 
lattice. Spatial dimension is normalized to lattice constant, distance between adjacent 
vortex centers (ref. 32). 
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Then equation (44) is substituted into equation (21) and the T f | A I'| 2 term 
is included. If the resulting solutions are forced to be periodic and the 
Gibbs free energy is minimized, the spatial variation of \xjj\ 2 can be 
calculated. The solution \\p\ 2 is found to have the symmetry of a square 
lattice in a plane perpendicular to the field direction (fig. 9). At the 
center of each vortex t \f vanishes and H = H a at these points. The lines 
correspond to constant \xjj\ 2 and constant H. Subsequently, Kleiner, 
Roth, and Autler (ref. 40) and Matricon (ref. 41) showed that a triangular 
lattice has a slightly lower energy (fig. 10). 



Gibbs Free Energy of Laminar Model and Vortex Model 

The ultimate theoretical test is to compare the Gibbs free energy for 
both models and determine which state is energetically more favorable. 
In each case an expression for the Gibbs free energy is calculated and 
then minimized. The H c \ can be determined by analyzing the dependence 
of G m i n on H. The minimum Gibbs free energy for the two models can 
then be compared using the H c i values. The pedagogical development 
that follows is largely due to DeGennes (ref. 38). 
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The analysis is begun by calculating a general expression for the total 
free energy of a superconductor. The Helmholtz free energy is the sum 
of the kinetic energy and the magnetic energy: 


^ s — E kinetic~^~ E magnetic"^- Eo 


The magnetic energy is 


^magnetic 


= ( hi 
J 8ir 


dr 


The kinetic energy is 


E k in e tic= j mv 2 (r)n s j dr 


The drift velocity v(r) is related to the current density by 

n s ev(r ) —j(r) 

The localized Maxwell equation is 


V x h= 



(45) 

(46) 


(47) 


(48) 

(49) 


where h is the local value of the magnetic field intensity. Then solving 
equation (49) for j(r), substituting into equation (48), solving for v(r) . 
and substituting into equation (47) yield 

Ekinetic J dr\]j'VX h] 2 (50) 

Finally, substituting equations (46) and (50) into equation (45) results in 

3^ s — E„ + ~ J (h z -h Ki\Vx h\ 2 )dr (51) 


In the following applications of equation (51), Eo — 0; that is, the normal 
state energy is taken as the zero of energy. 

Free energy of laminar model. — As illustrated by the laminar structure 
shown in figure 11, the local field h is determined by the London equation 
(eq. (6)). In one dimension with the assumed model, the London equation 
is 


388-980 0 - 70 -4 
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The solution in the superconducting region is 


h = H m 



cosh P 


( 52 ) 


(53) 


where P=d/2k and H m is the field in the normal region. Then when 
equation (51) is modified for the laminar geometry, the Helmholtz free 
energy is 



(54) 


Differentiating equation (53) with respect to x, squaring, and substituting 
the result and equation (53) into equation (54) yield 


r ‘ = £tii' d (i0[:J^ CO8 Ki0 + d^ 8tah '(iD] (55) 


When equation (55) is integrated, 



Figure 11. — Laminar model for mixed state. 
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1 H 2 m sinh 2 P 
SnP cosh 2 P 2 


(56) 


Using the multiple angle formula (sinh 2z = 2 sinh z cosh z) gives 


F s = 


H 2 m tanh P 
877 P 


(57) 


The free energy of the normal region is 


8 tt d 877 Pk 


(58) 


where #c=X/£. 

To convert the Helmholtz free energy to a Gibbs free energy, — 
must be added. Therefore, using the fact that B is the average of equation 
(53) in the superconducting region leads to 


BH^_ H flu re __ H I 2 

7T 477 477 \ d 


d/2 

dxH m 

0 


cosh ~ 
A 

cosh P 


HH m tanh P 
477 P 


Adding the three contributions to the Gibbs free energy of the lamina 
(eqs. (59), (58), and (57)) gives 

„ __H? n tanh P , HI 1 HH m tanh P 

laminar- 8v . p + 8 77 /V 4tT P (60) 


Minimizing G laminar with respect to H m yields H = Then if d is 4X/, 
or greater, tanh P ss 1 and 

G;nminnr = g ^(~ ■ff 2 + “^ £ ) (61) 

Now keeping in mind that the equilibrium state of a system corresponds 
to a minimum in the Gibbs free energy, Giaminar may be studied with 
respect to H. If H c < Vic, then G is a minimum when P is infinite. Infinite 
P implies infinite d or a complete Meissner effect. If H > H c IVk mini- 
mum G occurs for some finite P, which implies that it becomes energeti- 
cally favorable for magnetic field to be inside the superconductor. The 
laminar state becomes favorable when 

H — Hcii laminar) 

V K 


(62) 
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Free energy of vortex model. — The next step is to calculate H c \ for 
the vortex model. Once again, the Gibbs free energy is calculated and 
minimized to determine where the Meissner state breaks down. Near 
H c i the interaction energy between flux lines can be ignored. The Gibbs 
free energy of a material is equal to the number of vortices per square 
meter multiplied by the Gibbs free energy per meter of a single 

vortex minus the magnetic field energy density: 

Gvortex=B [f~ffj (63) 

where Z is the energy per unit length of a vortex. The energy Z can be 
calculated using equation (51): 

lfZ=\ drf |>2 + A2(VX fe)2] (64) 

J(r>€) OTT 

where If is the length of the flux line. This integration can be performed 
if the integral is first transformed to a surface integral (appendix C): 

l f Z=£ J da- [hx (Vxh)] (65) 

When integration is performed on equation (65) for the surface of a 
cylinder of radius the integrand is constant by symmetry and //cancels 
out: 


Z=£(27Tg)h(f)IVXh(f)\ (66) 


The h(i) 
equation 

where 


and |V X fe(f)| can be evaluated starting from the London 
h + \ 2 (V xVxh)= <p 0 d(r) (67) 


S(r) = 


0 

1 


r<£ 


If equation (67) is integrated over the area bounded by a circle C of radius 
r, which encloses the normal core, then 

J [h + \ 2 (VXVxh)]-da= j ^ O 82(r) , do - = <po 

By Stokes’ rule. 


( 68 ) 



MIXED STATE OF SUPERCONDUCTORS 


23 


J h • dcr + X 2 <J) |V X h\-dl = <p 0 (69) 

If k>r>%, then fh * da’—O because a negligible fraction r 2 /k 2 of the 
total flux is inside C. Then Vxfc may be assumed constant on C because 
of the cylindrical symmetry. Therefore, 


|Vxfc| = 


<Po 1 
27tX 2 r 


(70) 


The general solution of this equation is K 0 , the modified Bessel function. 
If V X h = — dh/dr , integrating by separation of variables yields 


= -^ln- 
2 tt\ 2 r 


(71) 


Now equation (66) can be evaluated at r = £ by substituting for h from 
equation (71) and for |V X h\ from equation (70): 


Z — 



(72) 


And now, finally, the Gibbs free energy of the vortex G vor tex (eq. (63)) 
can be calculated. If H < 4}7rZI<p 0 , the minimum in G VO rtex occurs at B — 0 
(type I superconductor). lf 'H>4irZl<p 0 i then the minimum occurs when 


dGyprtex 

dB 


■(£-£)-< 


and magnetic field penetration occurs when 


(73) 


H=H C1 = 


4-7tZ 

<po 


(74) 


Substituting equation (71) into equation (74) yields 

Hcl = 4S^ ln (i) (75) 

Because the original goal of this section is to compare the free energies 
of the laminar model and the vortex model, equation (75) must be alge- 
braically manipulated into a ratio of H c ijH c so that it can be compared 
with equation (62). This is done in appendix D and the result is 
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He l _ 7T In K 
He “ VM K 


( 76 ) 


Comparison of Gibbs free energy of the two models. — Now the Gibbs 
free energy of the two models can be compared. In figure 12, H c ilH c is 
plotted as a function of k for the functions derived in equations (62) 
and (76). The lower critical field H c \ is smaller for the vortex model than 



is 


o 



-2 1 I 1 i I 

0 12 3 4 


Ginzburg -Landau kappa, k 


Figure 12. — Curves H c \jH c against k for laminar model and vortex models are compared. 


for the laminar model. When H is greater than H c i (vortex) but less than 
H c i (laminar), the Gibbs free energy of the vortex state is less than the 
free energy of the Meissner state. And the Gibbs free energy of the 
laminar state is equal to the free energy of the Meissner state. Therefore, 
when H is between H c \ (vortex) and H ci (laminar), the Gibbs free energy 
of the vortex model is less than the free energy of the laminar model, 
and the vortex state is more favorable than the laminar state. DeGennes 
(ref. 42) has compared the states at higher magnetic fields and found 
that the vortex state is again more favorable. Goodman (ref. 43) has 
calculated the magnetization for both the vortex model and the laminar 
model. Figure 13 shows the comparison between these calculations and 
the magnetization measurements of Joiner and Blaugher (ref. 44) on 
Moo.ssReo. 15. Clearly, this result favors the vortex state. Subsequently, 
the vortex state has been substantiated by the NMR measurements of 
Pincus (ref. 45) and by the neutron diffraction experiments of Cribier 
(refs. 46 and 47). 
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FIGURE 13. — Goodman’s calculated magnetization curves for vortex and laminar states are 
compared with experimental data of Joiner and Blaugher (ref. 43.). Material, Moo. 85 Reo. 15 ; 
reduced temperature, t — TJT C — 0.522. 


The most convincing proof of the vortex state is the electron micro- 
scope pictures of U. Essmann and H. Trauble (ref. 48). They deposited 
iron on a lead-indium alloy in the mixed state. The flux lines in the alloy 
caused the iron to be deposited preferentially on the regions where the 
flux lines reside. An electron microscope picture of the result is shown 
in figure 14. A triangular lattice is formed by the flux lines in agreement 



FIGURE 14. —Electron microscope picture of triangular lattice of flux lines on surface of 
lead-indium alloy (ref. 47). Distance between flux lines, — 1000 A. 
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with the theoretical predictions. The flux quantum can be calculated by 
counting the dots and knowing the field of view and the magnetic field. 
The result is very close to 2 X 10~ 15 meter squared, tesla in agreement 
with predictions for the flux quantum. 

This experiment implies that it is energetically more favorable to 
create additional vortices as the magnetic field penetrates rather than to 
increase the flux in each vortex. This is theoretically agreeable because 
in equation (72) doubling <p 0 results in a line energy of 4Z compared to 
a line energy of 2Z if another vortex is added. It is interesting to note 
that the flux lattice has many properties of crystal lattices, plastic 
distortions, line dislocations, stacking faults, point defects, and holes 
(ref. 48). 


FLUX CREEP 

In the INTRODUCTION mention was made of the anomalous electric, 
magnetic, and caloric behavior of alloys that was observed by Mendels- 
sohn in 1934 and 1935. These experiments had the effect of focusing 
interest on alloy systems and lead to the discovery of the mixed state 
and type II superconductors. Some years later, Gorter (ref. 49) suggested 
that the Lorentz force resulting from electrical transport currents 
would cause regions of superconductivity to move and that this motion 
would result in energy dissipation. This concept along with the ideas 
of the mixed state and flux quantization set the stage for the discovery 
of flux creep and flux flow. 

Tube Magnetization Experiments 

Flux creep was first observed in tube magnetization experiments 
(ref. 50). The experimental setup, shown in figure 15, consists of a 
tubular sample with pickup coils inside and out. A probe for measuring 
magnetic field is in the bore of the tube. The magnetic field inside the 
tube H f is measured as an external magnetic field H is applied. Any 
flux that jumps into or out of the tube will be observed as transient 
induced voltages in the pickup coils. Abrupt changes in the flux in a 
material are referred to as flux jumps. Figure 15 also shows a plot of H ' 
as a function of H for a 3Nb-Zr alloy at 4.2 K. Induced currents in the 
sample shield the bore from the external magnetic field H up to 0.31 
tesla (curve A). Then the magnetic field begins to penetrate the bore 
and H t approaches asymtotically the value of H at very high field. 
When the external field is reduced (curve B), the magnetic field is trapped 
in the tube and the trapped field exceeds the applied field. The field H' 
becomes equal to zero at —0.31 tesla. When the same experiment was 
performed on a niobium powder sintered tube, the experimental data fell 
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External magnetic field, H, tesla 

Fjgube 15. —Experimental setup for tube magnetization experiments is shown and mag- 
netic field in the bore H of the tube is plotted against the field H outside the tube (ref. 49). 
Tube material, 3Nb-Zr; tube temperature, 4.2 K. 


on a straight line when the magnetic fields H and H ' were related by 


47 Twa c 

\M\ 


= B„ + H* 


(77) 


where w is the tube thickness, \Jl\ = // — #', and H*—i(H r +H). 
The B 0 and a c are constants* Kim, Hempstead, and Strnad were able 
to show that \M\ is proportional to the current in the sample, which 
leads to the critical state equation 

= J(B + B 0 ) (78) 

If J{B 4- B 0 ) > ot v , then magnetic flux enters the tube. 

Kim, Hempstead, and Strnad also made two other very interesting 
observations. If H is held constant, they found that \Jt\ decays loga- 
rithmically with time. Apparently the shielding currents decay in time, 
the magnetic field enters the sample, moves through it, and into the bore 
of the tube. They also noted that when the external field was being in- 
creased the flux entering the sample induced observable voltages in 
pickup coil Pi but not in P 2 (see fig. 15). When the magnetic field was 
being decreased, the flux leaving the sample induced observable voltages 
in pickup coil P 2 , but not in P u The smallest signals measured were 
estimated to correspond to 20 to 50 flux quanta. One possible explanation 


388-980 0 - 70 -5 
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for this observation is that flux enters the sample in unmeasurably small 
increments, perhaps individual flux quanta, but leaves the sample in 
increments of 20 to 50 quanta. 


Flux Motion Resistivity 

The observed decay of superconducting currents implies the existence 
of a resistance in the tube. The resistivity p of the tube can be derived 
by using a long solenoid approximation to get the current and by ap- 
proximating the tube with an equivalent LR circuit. The result is 


_ 2 rraw dH f 
P ~ H'-H dt 


(79) 


The existence of a mixed state resistance suggests that if a transport 
current is applied to a material in the mixed state, voltages should be 
observed. Kim, Hempstead, and Strnad (ref. 50) have in fact observed 
such voltages in NbZr. Figure 16 shows that these experimental data are 



Current density times magnetic field intensity, JH, (A/m 2 ) tesla 


Figure 16. — Voltage V is plotted as function of JH for curves corresponding to given cur- 
rent densities. Curve f is a plot of V against J(H-\-B 0 ) for all other curves (ref. 49). Alloy, 
3Nb-Zr wire; wire temperature, 4.2 K. 
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fitted by an equation of the form 

In V = J(H+ B 0 ) (80) 

Anderson and Kim (ref. 51) suggested a possible mechanism for the 
apparent movement of magnetic flux in type II superconductors. 
Anderson and Kim called the flux movement “flux creep” and postulated 
the following mechanism (see fig. 17). Flux quanta penetrate the sample 
individually when the sample goes into the mixed state (fig. 17(a)). But 
once inside, the flux quanta form into bundles (fig. 17(b)) which are ap- 
proximately the penetration depth in diameter. Inside each bundle the 
Abrikosov triangular flux lattice structure is postulated to exist (fig. 14). 
The postulate is consistent with the flux jumps observed in the pickup 
coils. They further postulate that the free energy of a superconductor 
varies spatially due to inhomogeneities, strains, dislocations, and other 
physical defects. These free energy variations provide pinning sites for 


(a) 






FIGURE 17. — Sketch of Kim* Anderson flux bundles and some other properties of flux creep 

mechanism (ref. 49). 
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the flux bundles. The strained regions have larger free energies than 
the unstrained regions (fig. 17(c)). These maxima in the free energy of 
the material represent minima for the free energy of the flux lines (fig. 
17(d)) because less energy is required to establish a normal core. The 
pinned flux lines therefore are surrounded by a free energy barrier. 
Anderson and Kim estimate the maximum free energy barrier to be 

A G max = (G n -G s )d* = (jjj?) £» (81) 

where the size of the inhomogeneity d is assumed to be about equal to 
£>. They further propose that a Lorentz type force acts on the flux bun- 
dles. This concept was suggested by the empirical critical state equation 
(eq. (78)). Flux creep occurs when the flux bundle is momentarily ther- 
mally excited to a free energy greater than the pinning energy. The 
Lorentz force then causes the bundle to move until it is caught by a 
pinning site (fig. 17(e)). Apparently, the thermal energy must then be 
lost to other flux bundles or to the crystal lattice; otherwise, the bundle 
would not become repinned. 

Anderson and Kim have also postulated that the rate R at which 
flux lines hop from one pinning site to the next is proportional to the 
negative exponential of the total barrier free energy. The total free 
energy is composed of two parts — the fraction p of equation (81) that is 
effective in pinning and the free energy associated with the Lorentz 
force. The Lorentz force per unit length of a flux line is 

f L = J -f (82) 

where c is the speed of light and <p 0 ~hc!2e. The force on a whole bundle 
is 


, .W J<p 0 n B V 
* L c c c 


(B3) 


where <p B is the flux in a bundle, n B the number of flux lines in a bundle, 
and U the effective length of a flux line (distance between pinning 
sites). The associated free energy Fl as a function of position x is 


F l 


JHkH'x 

c 


(84) 


This free energy of the flux bundle reduces the effective barrier free 
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energy that is pinning the bundle. The total effective barrier free energy 
Ft is then 

In this equation p is an adjustable parameter that ranges between zero 
and one to allow adjustment for the fraction of the maximum pinning 
energy that is effective. Next Anderson and Kim use the postulate that 
the number of flux bundles with a certain free energy decreases expo- 
nentially with the free energy. This implies that the thermal activation 
of flux lines in the bundle over a free energy barrier is governed by a hop 
rate R of the form 


R ~ 0) 0 €~ F T^' t 


( 86 ) 


In this equation a> 0 , which is the vibrational frequency of the bundle, is 
estimated to be between 10 5 and 10 10 hertz. The flux bundles move in 
the direction of the gradient of magnetic pressure a, which may be 
defined as 


a— V 



jxH 

c 


(87) 


Then if the pinned flux bundles are viewed as a supply of flux lines, 
the flux entering and leaving, a small volume obeys the diffusion equation 



These equations ((85), (86), (87), and (88)) can be combined to study flux 
motion in the tubular geometry of the tube magnetization experiments. 
In one dimension the diffusion equation becomes 


dH __ d ( <poR \ 

dt~ dx \ k ) 


and the magnetic pressure is 


FLdH 

47 T dx 


(89) 


(90) 


Then taking the time derivative of a and substituting for dH/dt from 
equation (89) yield 
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da = 1 (podfjdR | H <podVt 
dt 4 itt X dx dx 4 >it X dx 2 


If equation (86) for R is written in terms of a (by substituting eq. (87) 
into eq. (85) and then substituting the result into eq. (86)) and the result 
is substituted into equation (91), then 


da _ ^ dH de a ' a 1 , ^ „ d 2 e a/a , 
dt dx dx dx 2 


where 


If the magnetic field is assumed to be uniform, then dHjdx = 0 and the 
first term on the right side of equation (92) vanishes. The remaining 
differential equation has a particular solution of the form 


a = ai (in 


aix 2 + bx + c 
2 K 0 


— In t'j 


This result agrees with the logarithmic time decay of the magnetic field 
trapped inside tubes of Kim’s tube magnetization experiments (ref. 52). 
Furthermore, if equation (86) is solved for a (as described in the derivative 
of eq. (92)) and it is assumed that there is some critical R c which cor- 
responds to the lowest hop rate that can be observed experimentally, 
then an equation very much like equation (78) can be derived. 

Kim and Anderson also derived the experimental result equation (80) 
and the power dissipation from Faraday’s law: 


V x £ = 


In one dimension. 


z (po dR. <p 0 a)o d 

ck dx kc dx 


( e -F 0 lkT e al< 


where F 0 = p(//?/8tt)^. Integrating with respect to x gives 


(96) 
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The power dissipated P is then 


(97) 

The existence of power dissipation is consistent with the severe thermal 
instabilities observed in superconductors that carry high current density 
in large magnetic fields. 

Kim, Hempstead, and Strnad (ref. 49) have found that in the limit of 
large values of a the exponential hop rate law (eq. (86)) is no longer 
applicable. Instead, R is proportional to a. This might be expected be- 
cause for large a, F 0 — (ajai)kT becomes comparable to kT. This new 
regime is called the flux flow state. 

FLUX FLOW 

Early Flux Flow Experiments 

Kim, Hempstead, and Strnad first observed the flux flow state as a 
deviation from the flux creep results (ref. 49). In figure 16 the composite 
curve f would be a straight line if flux creep were the only cause of the 
voltage. But at larger values of a, the curve bends over indicating that 
a new mechanism is contributing to the voltage. They measured the 
voltage-current -magnetic-field characteristics of several different 
materials. Figure 18 shows some examples of their results for NbTa and 
Pbln. These measurements were made with the magnetic field perpen- 



(a) Two samples of Nb 0 .s 
Tao.s with different 
amounts of pinning. 
Temperature, 3 K; mag- 
netic field, 0.2 tesla. 


(b) Comparison for three dif- 
ferent magnetic fields and 
also for plated and un- 
plated samples of Pbo .83 
Ino. 17 . Temperature, 2 K. 


FIGURE 18.— Flux flow voltage is shown as function of current (ref. 53). 
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dicular to both the surface of the samples and the applied current. Volt- 
age appears when the current is greater than some magnetic-field- 
dependent onset values. The voltage is directly proportional to the 
current above the onset current. The voltage is also observed to be an 
increasing function of the magnetic field (fig. 18(b)). At temperatures 
sufficiently below the critical temperature T c , the ratio of the flux flow 
resistivity p/ to the normal resistivity p n obeys the simple relation 


pl=JL 

pn H C 2 


(98) 


This relation is very suggestive of a law of corresponding states. A re- 
sistivity like pf would be observed if the current flowed uniformly through 
the superconducting regions and the normal cores of the flux vortices. 

At first thought it seems unlikely that the current could flow uniformly 
through both the superconducting regions and the normal regions. The 
superconducting regions should simply electrically short the normal 
regions. If, however, the flux flow state is similar to the flux creep state, 
then the flux lines move under a Lorentz force. Also, the motion might 
generate electric fields that drive current through the normal cores. This 
would result in a resistivity that is proportional to the number of flux lines 
inside the superconductor. The model just described was first proposed 
by Kim, Hempstead, and Strnad (ref. 53). 

The flux flow voltage can be expressed as a function of J and H . Figure 
19 shows a type II superconductor with the magnetic field perpendicular 
to the current and to the surface of the superconductor. The vortices are 
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Figure 19.— Flux flow state. 
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indicated by simple current loops. The magnetic flux density equals the 
density of flux lines multiplied by the flux quantum: 

B = mpo — H (99) 

The magnetic intensity and the current density are related by Ampere’s 
law: 

vxH=i^-=i<p 0 ^=— J (ioo) 

dy r dy c 


The magnetic pressure is 

A /£ 2 LM = I /4r ZjXHJ^meo 

dy \87 t) 47 t dy 47 t\c j ) c c 

And if the magnetic pressure is distributed evenly over 
the Lorentz force on each flux line per unit length is 

f L = J -f (102) 

When the Lorentz force is greater than the pinning force, the flux line 
moves. It is further postulated that the flux lines experience a viscous 
drag force. The flux line accelerates until the viscous drag force balances 
the Lorentz force: 

fL^y]VL (103) 

where t) is a viscosity parameter and Vl is the terminal velocity of the 
flux line. The power P dissipated per unit volume by the viscous flow is 

P = nf L v L = - (104) 

7} C 

Therefore, it is possible to derive the voltage per unit length V from 
equation (104) for P — VJ and H = mp 0 : 


J (101) 

all the flux lines. 


V=~JH (105) 

7)C 2 

The linear dependence of V on J is in agreement with the experimental 
results shown earlier (fig. 18). The magnetic field and the temperature 
dependence is shown by figure 20. The flux flow resistivity p/= VjJ is 
directly proportional to the magnetic field when T < T c . 
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FIGURE 20, — Normalized flux flow resistivity plotted against magnetic field at several 
reduced temperatures for Nb 0 .i Tao.9. Resistivity also shown plotted against temperature 
at several values of magnetic field. Flux flow resistivty is normalized to normal state 
resistivity (ref. 54). 


Other Flux Motion Experiments 

There are other experiments which support the concept of flux 
motion and the flux flow voltage. Pearl (ref. 55) proposed to test the 
concept of a potential difference resulting from flux motion. He proposed 
to produce a continuous motion of flux lines by some external means 
and to look for a voltage. Later (ref. 56) he built such an apparatus 
(fig. 21). Magnetic field lines extending from the screw to the iron shell 
could be moved through the superconductor by rotating the screw. No 
transport current was passed through the sample, and yet voltages 
were observed across the sample that were roughly proportional to the 
rate of rotation. 



Iron shell 
Magnetized screw 
^-Superconductive 


Direct current 
voltmeter 


Mixed state 
area 


Motion of 
magnet 


Figure 21. —Apparatus used by Pearl to investigate flux flow voltage (ref. 56). 
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Hans Meissner (ref. 57) performed another experiment which tends 
to confirm the concept of the motion of flux lines under a transport cur- 
rent induced Lorentz force. Meissner introduced a flux bundle on one side 
of a tin (Sn) film and observed it some time later with a pickup coil 
on the other side of the film. The flux flow velocity is determined by 
measuring the transit time and the film width. Meissner has observed 
flux flow velocities between 3 X 10 2 and 50 X 10 2 meters per second 
with current densities up to 10 9 amperes per square meter. 

Another experiment on thin films by Giaever (refs. 58 and 59) con- 
tributed results which were consistent with the flux flow model. Giaever 
investigated the properties of two parallel thin films, 1000 A thickness, 
in a perpendicular magnetic field. When he passed current through the 
primary film, he observed voltages across both the films. Figure 22 
shows the model he proposed to explain these results. Both films are 



Figure 22.— Flux sharing model proposed by Giaever to explain his experimental results 

(ref. 59). 

assumed to share the same flux bundles. When the transport currents 
in the primary film are sufficient to overcome pinning, the flux moves in 
both films. Therefore flux flow voltages are observed in both films. 
Some of Giaever’ s results are shown in figure 23, which shows the 
voltage in the primary film when the current is applied to both the 
primary and the secondary films. The family of curves shows that the 
current in the secondary apparently exerts a force on the flux line which 
can either add to or subtract from the force exerted by the current in 
the primary. 
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FIGURE 23. — Voltage across primary film plotted against primary current for various cur- 
rents in secondary film (ref. 59). Temperature, 7 1 — 3.65 K. 


Theoretical Basis for the Flux Flow Voltage 

The exact nature of the source of the flux flow voltage has caused consid- 
erable controversy. The fact that p/ is proportional to H , or the number 
of normal cores (eq. (98)), suggests that the voltage is caused by transport 
current flowing through the normal cores of moving flux lines. However, 
Hempstead and Kim (ref. 60) suggested that the flux flow voltage was an 
induction voltage. Several authors (refs. 61 and 62) have argued that in 
the standard flux flow experiments the flux lines move across the sample 
but that dB/dt = 0 and therefore no induction voltage can exist. Other 
authors (refs. 63 and 64) argued that dB/dt = 0 means only that V X E = 0 
and not that E = 0. 

Goodman (ref. 65) has successfully derived the flux flow voltage from a 
local version of the second Ginzburg-Landau equation. Ignoring the 
first term on the right side of equation (22), the starting equation is 


V 2 a = — 122 
c 



a . 

-2e™) 


(106) 


where a is the local value of the magnetic vector potential, j is the local 
current density, i{f=i[f o e i0 , and VO is defined by the relation 
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dl * V0 = Ittub 


( 107 ) 


where n B is the number of flux quanta per bundle. He also defines the 
average density of vortex lines N and an incremental area dS. It is as- 
sumed that dS is much larger than the distance between flux lines. Then 
taking a line integral of the Ginzburg-Landau equation yields 

<f j ■ dl = — (j a-dl-f ~-V6-dl) (108) 

C Ids me - \ Jds Ids 2e ) 

Converting equation (108) to surface integrals using Stokes’ theorem gives 


4 TT 

C 


f (Vxj) 

JdS 


dS = 


1 6 ye z tyl ~ f 
me 2 LJds 


(V X a) • dS 


iV«Po] 


(109) 


Recalling that xjj'i = ns and = mc 2 l^7rn s e 2 and substituting yield 


H + — (Vxj)= iVijOo 

c 


( 110 ) 


where H and J are the mean field and current density. Equation (110) 
can be written in either of the following forms using Ampere’s law 
(Vxff = (4ir/c)J): 


H- kl\7 2 H^N<p 0 (111a) 

or 

/-XfVU=g(VXiV) (111b) 

Thus, if N is known, J and H are specified. Next Goodman assumes 

a vortex line continuity equation: 

Vx(u,XJV)=^ (112) 

at 


Then taking the time derivative of equation (110) gives 


dH 

dt 


+ 4>irk 2 




(113) 


Substituting equation (112) and Faraday’s law of induction (Vx£ = 
— (1 /c) (dH/dt)) into equation (113) yields 
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-V X E +^|[v x =&- [V X (»,. X TV) ] (114) 

Equation (115) is obtained by eliminating the curl on both sides of the 
equation and setting djjdt — 0: 


E=—Nx vl) (115) 

c 

Goodman has ignored any possible constant terms and gradient terms 
in equation (115). Equation (115) is identical to the electric field that is 
measured by a stationary observer in a coordinate system moving with a 
changing magnetic field (ref. 66). 

THEORIES OF THE MOTION OF VORTICES IN SUPERCONDUCTORS 

Before going into the two major theoretical models that have been 
proposed for flux flow, it is possible to gain some insight into the problem 
by a simple analysis of the change in momentum resulting from the 
motion of a single vortex. Figure 24 is a sketch of a normal core (radius 
r) moving with velocity v L . Outside the core the superconducting elec- 
trons of the vortex are assumed to have momentum: 



When the vortex moves as shown, the superconducting electrons on 
the front boundary of the core are converted to stationary normal 
electrons. The resulting momentum loss in the boundary layer is ap- 
proximated by 




FIGURE 24. — Normal core (radius r) moving with velocity \ l . 
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dp __ (Number of electrons)(Momentam change) ( Trrn s ds)pe 
dt Time interval dt 


(117) 


Setting VL — dsjdt and substituting equation (116) into equation (117) 
gives 

^% = ~\ nn ^ lVL (H8) 


Finally, substituting for h from <p 0 = hcj2e yields 


dp _ 1 n s e 


dt 


2 c 


Vl<Po 


(119) 


Because an equal amount of momentum must be added on the right side 
of the vortex, the resulting momentum lost must be twice equation (119). 
The theories of Bardeen-Stephen and Nozieres-Vinen (which follow) 
postulate mechanisms for this momentum loss. The rate of change of 
momentum corresponds to a Magnus type force (ref. 67), which acts on 
the vortex and is perpendicular to its direction of motion. 


Bardeen-Stephen Theory 

Two theoretical models have been proposed which attempt to explain 
the phenomena of vortex motion in type II superconductors. The theory 
of Bardeen and Stephen (ref. 68) finds that, if the transport current 
causes a force on the flux line, the resulting motion generates an electric 
field that drives the current through the normal core of the vortex. 

Bardeen and Stephen assume that the Lorentz force J<p 0 lc is the 
driving force on the flux lines. They assume that the most important 
scattering process is electron-lattice scattering rather than electron- 
electron scattering. An isolated vortex is considered; the temperature 
is near 7=0 K and the mean free path is assumed to be greater than 
the coherence distance. Bardeen and Stephen assume that the flux line 
contains a normal core. This means that the energy gap A is equal to 
zero inside the core. Outside the core there is assumed to be a tran- 
sition region where the energy gap increases from zero to its equilibrium 
value. They further assumed that the transport current density and the 
normal current density generated by the motion are small compared 
with the supercurrent flow just outside the core. Then they studied the 
current distribution in the transition region starting with the equation of 
motion from the local London theory: 


dv s 

dt 


„ e 

— — V/x 0 H — 

^ m 



( 120 ) 
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where fi 0 is the chemical potential per unit mass with zero current and 
field. They succeed in deriving expressions for the components of the 
velocity of the normal and superconducting electrons. Combining these 
equations with calculations of the power dissipation inside and outside 
the core, they obtain several interesting results. First, the velocity of 
electrons inside the core is equal to the velocity of the electrons in the 
transport current. If the core radius r 0 is assumed to be related to H c 2 by 


f P° —i j 

_ _ — n C 'i 


H 


2 irrl 


( 121 ) 


the law of corresponding states (eq. (98)) can be derived. They also con- 
clude that the Hall angle 6h is equal to the Hall angle of the normal ma- 
terial at the magnetic field in the normal core; therefore, 

tan 6h = (—)h ( 122 ) 

\mcj 

where r is the normal electron relaxation time. This same result was 
obtained earlier by Miller (ref. 69) using a two fluid model. 

Nozieres-Vinert Model 

The other theoretical model was proposed by Nozieres and Vinen 
(ref. 70). Nozieres and Vinen proposed that the motion of flux lines is 
analogous to the motion of vortices in liquid helium. The vortex of super- 
conducting electrons is subject to a Magnus force, which is perpendicular 
to the direction of vortex motion. DeGennes and Matricon (ref. 71) had 
proposed a similar concept earlier. They considered the force on a flux 
line with an applied transport current. In a reference frame moving at a 
constant velocity Vi ,, superfluid electrons represent the current 

j~n s {-e)(-v L ) (123) 

Therefore, the force per unit length on the flux line is 

f l =j^=-n s e(^j <Po ( 124 ) 

The Nozieres and Vinen model incorporates the following assumptions. 
The penetration depth is much greater than the coherence length. The 
material is homogeneous, with the mean free path much greater than the 
coherence length. The temperature is near 0 K so that no normal excita- 
tions exist outside the core. The velocity distribution of the other vortices 
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is assumed to be constant over the core of the vortex, with a sharp bound- 
ary assumed for the core. Far outside the core the London equation is 
assumed to be obeyed. The entire development is carried out in the 
crystal lattice frame of reference to avoid problems with the effect of the 
positive background. Then, using the London and Maxwell equations 
together with Euler’s equation for a non viscous charged fluid (ref. 33), 

v) *’' ir + 2 ~ n) 

— {v s — vi) X (V X v s ) (125) 

where v s is the total superfluid velocity and vl is the velocity of the flux 
line, they derive Bernoulli’s theorem for the charged fluid: 

/x + eF+-^ (v s “ Vl) 2 = Constant (126) 

Next they derived the force on a line by considering the flow pattern in 
the vortex frame. To do this they constructed a cylinder of radius r > g 
which was concentric to the flux line. Then, summing the electromagnetic 
forces, the fluid pressure, and the forces caused by the net momentum 
flowing through the surface of the cylinder yielded the total force as 

Ne 

F=~ {vsi~-v L )<p 0 (127) 

where v s i is the local superfluid velocity in which the line is moving. 
This force per unit length is the Magnus force. 

The concept of using the equations of fluid mechanics to describe 
superconductivity is not new, it was first proposed in 1933 (ref. 17). 
London derives equations (125) and (126) in a discussion of these earlier 
attempts (ref. 33). 

Then when the vortex core is approximated by a normal cylinder with 
radius a, the electrons in the core are assumed to have a uniform drift 
velocity i ? nc . At the boundary between the core and the remainder of the 
vortex there is a discontinuity in the superconducting wave function. 
Following Bardeen and Stephens, Nozieres and Vinen assume v n c = v s \ . 
They also assume that the Magnus force is divided evenly between the 
region outside the core and bulk of the core. The surface part of the 
total force would come from an electrostatic contact potential at the inter- 
face. They then assume that no contact potential exists and set the re- 
maining half of the Magnus force equal to the rate that momentum is 
being lost in the core. This gives the equation of motion of the normal 
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cores: 

l — (»,1 - V,) X ip—Nm — 7 TO 2 = 0 (128) 

l C T 

This equation may be compared to the Bardeen-Stephen result if the 
relation between core radius and field is 


2ira 2 


H ( 


c2 


(129) 


When this is substituted into equation (128) the resulting equation is then 
solved for the resistivity and yields the law of corresponding states 
(eq. (98)). The Hall angle that results from the Nozieres-Vinen analysis 
is independent of magnetic field: 


tan 6 h = — r (130) 

me 

Comparison Between Theoretical Predictions and Experimental Results 

An obvious experimental test of these two theories is to measure the 
Hall angle and determine whether or not it has a magnetic field depend- 
ence. Figure 25 shows the Hall angle predictions of the two models. The 

Niessen and 



Figure 25. — Theoretical predictions of Bardeen-Stephen and Nozieres-Vinen are compared 

with experimental results. 
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presently available experimental data are also included in the figure. The 
results of Niessen and Staas (ref. 72) on dirty (mean free path < £ 0 ) 
NbTa does not agree with either theory. Very recently a similar experi- 
mental result was obtained by Hake (ref. 73) on TiMo. Hake also finds 
the Hall angle in the mixed state to be greater than the Hall angle in the 
normal state. However, because both theories assume pure materials 
these results are interesting but not surprising. The Hall angle results of 
Reed, Fawcett, and Kim (ref. 74) on high purity niobium (resistance 
ratio, R300/R4.2 “ 1550) agree more closely with the Bardeen-Stephen pre- 
diction. However, other experiments (by Fiory and Serin) on niobium 
(resistance ratio, #300/^4.2— 3000, refs. 75 and 76) pointed out the 
effects of pinning. Both theories assume that there is no pinning. There- 
fore, in order to distinguish between the two theories the effects of pin- 
ning must be eliminated. Fiory and Serin first proposed to do this by 
de-pinning the flux with a 20-kilohertz, 1.50 X10” 2 tesla magnetic field 
superimposed on a constant magnetic field. The direct current density 
used was 1.350 X 10 7 amperes per square meter. They obtained results 
(ref. 75) which more closely agreed with the Bardeen-Stephens theory. 
In their second paper (ref. 76) they used direct current densities of 
1.1 X 10 8 amperes per square meter to overcome the pinning. Their re- 
sults this time agreed with the Nozieres-Vinen prediction. By way of 
comparison, the maximum current density used by Reed, Fawcett, and 
Kim was 1.090X10 7 amperes per square meter. Still another experi- 
ment performed on niobium by B. W. Maxwell (ref. 77) supports the 
Nozieres-Vinen model. Maxwell measured helicon resonance parameters 
in the mixed state. Relating these measurements to the Hall angle he 
found it to be field independent. The present experimental evidence 
seems to support the Nozieres-Vinen model when pinning effects are 
overcome, whereas the Bardeen-Stephens model applies when pinning 
effects are still effective. Nozieres and Vinen (ref. 70) suggest that an- 
other crucial experiment would be the measurement of the dispersion of 
helicon waves in very pure materials (o) c r> 1). The theoretical predic- 
tions are 


o> = 


FA 2 


PA 2 


eH 

me 

eHJL 

me H C 2 


Nozieres-Vinen model 
Bardeen-Stephens model 


(131) 


where <0 is the frequency of the helicon waves. This experiment must 
await higher purities than are presently available. 
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RIGIDITY OF FLUX LATTICE 

The Abrikosov flux structure has a rigidity that results from the 
repulsive force between flux lines. DeGennes (ref. 42) has derived the 
force between two vortex lines. He starts with the assumption of a 
modified London equation: 

X 2 (V x v X h) = (p 0 [8(r—ri) + S(r-r 2 )] (132) 

The delta function terms are added to include the contributions to the 
field resulting from two flux lines at r\ and r 2 : 

h=h 1 (r) + h 2 (r) (133) 

The free energy of this system is given by equation (51), and equation 
(65) is 

F=^ j hx (V x h) da 
X 2 r 

= ftoj (dw + dot) * (hi + h 2 ) X (Vxfcj-f Vxfe 2 ) (134) 

It is understood that the surface of integration is the cylindrical surface 
of the two cores (|r—r/|=f). The indicated vector operations in the 
integrand of equation (134) result in eight terms: 


^ [ JdcTi • hi X (V X jfej) + * hi X (VXfe 2 ) + J do-i * fe 2 x (Vxfej) 

+ fd&i 1%2 x (Vxfc 2 ) ~f fd<r 2 • hi X (V X hi) 
+ Jdcr 2 -/ii X (V X fc 2 ) + J d<r 2 ‘ h 2 X (V X h x ) 
jd&2 * h 2 X (V X fi 2 ) ] (135) 

The first term and the last term represent the self-energy (eq. (65)) of 
each vortex fine. DeGennes assumes that the current (V X h ) at one 
vortex core due to the other vortex can be neglected. Therefore, in equa- 
tion (135) the second, fourth, fifth, and seventh terms are dropped. The 
remaining two terms represent the interaction potential between the two 
flux fines t/i 2 : 
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From equation (70) and its general solution K 0 {x) the modified Bessel 
function of order zero, 


Vx& 2 


<Po 1 


and 


27r\ 2 I r-~r 2 I 

Therefore, 

Jh.x(V X h 2 ) • d<r 2 = j ^K, 

Jh, X (V x M • rf«r 2 = 2 tt (^) 2 K 0 ( |ri ~ r ^ 1 ) 

<Po „ / |ri-r»l\ 

2ir\ 4 Ao V X / 


(137) 

(138) 


(139) 


(140) 


The other term in equation (136) gives an identical result. Therefore, the 
interaction energy is 

(14I) 

Substituting the field at vortex line 2 due to vortex line 1 (eq. (138)) into 
equation (141) results in 

Ui2=^h l2 (142) 


The modified Bessel function 
for large ri 2 and 

for small ri 2 . Figure 26 shows a plot of X 0 (ri 2 /X): note that it becomes 
very small at something less than ri 2 /\ = 3. The force exerted by one 
flux line on the other can be calculated from the potential U\%\ 
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X 

FIGURE 26.— Modified Bessel function plotted against its argument (ref. 78). 


fix 


dUn 

dx 


<po dh\2 

47 T dx 


(143) 


where x = n 2 /X is taken as one axis in a rectangular coordinate system. 
DeGennes then identifies dh^ldx as — (47 r/c)y v , which is just the current 
density predicted by Maxwell’s equation V X H= (torlc)J. This gives 
a force 


fix^irjy (144) 

which is a Lorentz type force. The flux lines are therefore forced apart, 
but they cannot leave the superconductor because if they did more 
energy would be required. The additional energy would be needed to 
create surface currents large enough to exclude more magnetic flux. 
The repulsive force between flux lines therefore results in rigidity of the 
flux lattice. Kim, Hempstead, and Strnad’s theory (refs. 49 and 53) of flux 
flow implicitly assumes that the entire flux lattice moves as a unit. 

DeGennes’ result (eq. (144)) is essentially valid when the flux lines are 
more than a penetration depth apart. If the flux lines are less than a 
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penetration depth apart, then all the terms of equation (135) which were 
discarded must be included. 


VISCOUS DRAG FORCE 

In the section on FLUX FLOW, a viscous drag force on the moving 
flux lattice was postulated (p. 33). It was needed to explain why the flux 
lattice does not accelerate indefinitely. The proposed mechanisms for 
the viscous drag force are discussed in this section. In actuality some 
combination of these effects, perhaps even all of them, may contribute 
to the viscous drag force 

/= yvi, (145) 

Before proceeding with the various models for the viscous drag force, 
it is useful to derive an empirical expression for the viscosity parameter 7 ). 
This result can be compared with the results of the various theories to 
evaluate them. Dividing equation (105) by J yields the flux flow resistivity 


Pf 




7] C* 


(146) 


Then, substituting this for p/ in the law of corresponding states (eq. (98)) 
and solving for r\ yield 


V empirical 


(ppH C 2 ( 0 ) 

pnC 2 


(147) 


Equation (147) can be written in several alternate forms (ref. 53) using 
the expression for H c 2 ( 0 ) derived by Maki (refs. 79 and 80) and using 
several relationships between the fundamental parameters: 


or 


7 ? empirical — TthN (0) A# 

= 6.4 X 10~ 12 yT c 


(148a) 

(148b) 


3 <T n <pl 

27 T 2 C 2 £ 2 


(148c) 


where 7V(0) is the normal state density of electron states, y is the elec- 
tronic specific heat constant, and cr„ is the normal conductivity. Equation 
(148) points out that 7] empirical depends only on N(0) and A 0 . Notably it is 
independent of the mean free path of the electrons. 
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The first viscosity mechanism to be discussed is normal eddy-current 
damping, which was first proposed by Volger, Staas, and Van Vijfeijken 
(ref. 81). If equation (145) is substituted into equation (104) for F/ v , the 
power dissipated is 


P = n7)vl (149) 

At low temperatures where the only normal electrons are in the core, the 
total power dissipated should be equal to the number of cores n times the 
power dissipated in each core: 

P^n$SaE*dxdy (150) 

As mentioned previously, when flux moves there is an associated elec- 

tric field (eq. (115)). When the applied current and field are perpendicular, 
E-— ( viJc)H . An expression for r) is obtained by substituting this into 
equation (150) on the right, substituting equation (149) into equation 
(150) on the left, and solving for tj: 

Yi= ^SI H ^ dxdr - E ^ <i5i » 

Maki’s expression (refs. 70 and 71) for H c 2 (0) can be manipulated into 
the form 


#c 2 (o) =2™.^ ( 152 ) 

Then if equation (151) is multiplied and divided by H* 2 

h&sH' < 153 > 

Kim, Hempstead, and Strnad (ref. 53) have pointed out that equation 
(153) is only equal to equation (147) near H c2 , and the 7} given by the 
normal eddy-current damping model is always much less than the 
empirical value. Therefore, this mechanism may well exist, but it is by 
no means the primary source of viscosity. 

Stephen and Bardeen (ref. 82) proposed a cycloidal damping scheme. 
They pointed out that if a flux line does not move too fast,nhe entire 
vortex structure moves as a unit. Under these conditions the electrons 
in the vortex move in a circular path when the flux line is stationary 
u$(r). But when the vortex is moving, the electrons move on a cycloidal 



MIXED STATE OF SUPERCONDUCTORS 


51 


path v s {r~ VLt). This situation is somewhat analogous to the motion of 
an electron in crossed electric and magnetic fields. Stephen and Bardeen 
argue that likewise in the case of vortex motion an electric field must be 
present to exert the force that results in the cycloidal motion. If 
R = ( r — VLt ) , the acceleration experienced by the vortex electrons is 

(154, 

If the flux flow motion is steady state, dvijdt — 0. Because the vortex is 
assumed to be rigid, dr/d£=0. The electric field associated with the 
force on the vortex electrons is 

(£)E=(v l -\7)v s (155) 

A detailed analysis of this electric field in the region of the vortex core 
indicates that inside the core the field is constant and equal to 


(156 > 

where p is the momentum of the paired electrons defined by the quantum 
condition 


2 



dl = h 


Then 


P 


A 


(157) 


(158) 


and £ is the radius of the circular path of integration. Differentiating 
equation (158) and substituting into equation (156) gives 


E c =— 


vl h 
e 2? 


(159) 


Now a viscosity parameter can be calculated by equating the total rate 
of dissipation to the dissipation in the core: 


rjvl — iT^a-niEl) 


(160) 
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Then if equation (159) is substituted into equation (160), solving for V 
yields 


_ rrcrnh 2 

V ~ eHe 


(161) 


Multiplying and dividing by <p% = fe 2 c 2 /4e 2 in equation (161) gives 


<Tn<fl 

71 4nc 2 g 2 


(162) 


Equation (162) is 7 t/ 6, or roughly one-half the empirical viscosity (eq. 
(148c)). 

Tinkham (ref. 83) has suggested an additional source of damping. It 
results from the variation of ip when the flux lines move. At a point inside 
a material which is in the flux flow state, the superconducting wave 
function oscillates between \p 0 (superconducting state) and 0 (normal 
state) as normal cores pass through the point. Tinkham approximates 
the frequency of this oscillation by (oszvijg. The superconducting elec- 
trons thus gain and lose the condensation energy w/27t times each second. 
If the instantaneous value of ip lags behind the value appropriate for the 
moving Abrikosov structure by a time r, then a fraction cor of the energy 
is dissipated. This energy goes to the production and/or heating of normal 
electrons. Tinkham approximates the Gibbs free energy by the second 
term of equation (33), which is the Ginzburg- Landau expression. Then 
equation (163) is obtained by multiplying and dividing the approximation 
by ipg (eq. (36)): 

r = Hj{T) ip 0 (H, T) 

47 T xp 0 ( 0 , T) 

For convenience, 


(163) 




MH, T) 2 
MO, T) 


The power dissipated per unit length is the amount of energy lost per 
cycle times the number of cycles: 


P = (Go>t) o>it £ 2 = X 2 j p 2 tt£ 2 (164) 

Again the plan is to set the power dissipated to r\vl and solve for 77. But 
before doing this Tinkham approximates x 2 and r. He assumes that the 
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relaxation time is limited by the uncertainty principle 

ft h (1 + f 2 ) 

T A A 0 (1-f 2 ) 


(165) 


where t = T/T c or the reduced temperature. He further approximates 
X 2 by 

X 2 = 1 — 6 2 (166) 


wher eb = B!H C 2 - Then substituting equations (165) and (166) into equation 
(164), setting P equal to i)vh r} can be solved for 


7, = T f X 2 W,T) 


nm 

4A„ 


(1-i 2 ) 


■(.I t l 2 ) 

(1 ~* 2 ) 


(167) 


This expression is quite different from riempincai • It has both a tem- 
perature and magnetic field dependency. However, the temperature 
and magnetic field dependence in r)empirtcai was ignored. Therefore, it 
is possible that equation (167) is, in fact, more exact. If equation (167) is 
substituted into equation (146), 


Pf 


4 Hyp 
c 2 tH 2 X 2 


(168) 


Then for dirty (mean free path £ 0 ) superconductors 

1 ^c 2 h 
pn 47 T A 0 A 2 (0) 

U rHhX 2 (T) 

C2 <Po 


and 


Substituting these into equation (168) gives 

Pf 4 H h XHT ) 1 

Pn ir. H ci {T) A 0 t(T) \ 2 (0) xHH, T) 


If \ 2 (0 )/\ 2 (T) is set equal to (1— t 4 ) and equations (165) and (166) are 
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substituted into (169), 


= 4 1 b 

Pn 7T (1 + fi) (1— i 2 ) 


( 170 ) 


Tinkham (ref. 83) then adds this to his own generalized version of the 
normal eddy current damping model and obtains 


p n (l + t*)(l-b*)+b*+(b-b*)t K J 

Note that at £ = 0, pflp n ~HIH C 2 and at H~H C 2 , P//Pn = 1. This relation 
can be compared directly with the experimental results (see fig. 27). 
The agreement is remarkably good, especially at low temperatures. 



FIGURE 27. — Flux flow resistivity normalized to normal resistivity is plotted against applied 
magnetic field normalized to upper critical field H c z. Tinkham’s (ref. 83) results are 
compared with experimental results. (Temperature in Kelvin degrees represented by 

T.) 


It has been suggested by Kim, Hempstead, and Strnad that the mecha- 
nisms of Bardeen-Stephen and Tinkham are independent (ref. 53). Their 
relationship may be better understood when experiments are performed 
in sufficiently clean materials. The Bardeen-Stephen model should be 
dependent on the mean free path, because it depends on collisions in 
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the vortex core. The Tinkham result is independent of mean free path. 
Both predictions are the correct order of magnitude to explain the 
experimental results. 


FLUX PINNING MECHANISMS 

In the FLUX FLOW section, the onset of flux flow was described. 
When the Lorentz force is sufficiently large the flux lattice moves and a 
voltage is observed. It is thought that a minimum Lorentz force is nec- 
essary to overcome the pinning forces exerted by crystal imperfections 
on the flux lines. The nature of these pinning forces is discussed in this 
section. The primary sources of pinning are thought to be dislocations 
and impurity atoms. Figure 28 attempts to show schematically several 
different types of pinning mechanisms. In sintered materials such as 



Superconducting 

material 


Figure 28 . - Cross-sectional view of three different pinning mechanisms. 


Nb 3 Sn, where the sintering is incomplete, precipitates and cavities are 
thought to be the most important pinning mechanism (ref. 84). A flux 
line will become entwined in the defects (fig. 28) such as to minimize 
its energy taking account of the distance between flux lines and the size 
and density of the impurity atoms and cavities, Friedel, DeGennes, and 
Matricon (ref. 85) considered the case where one line is pinned by each 
cavity. They estimate the pinning force F p to be 


F p = 



(172) 


by equating F p to the line energy (eq. (72)) and setting k~d\ the distance 
between two pinning sites. If the density of cavities and/or precipitates 
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is such that it becomes likely that each flux line is pinned by more than 
one cavity, then presumably F p should be divided by the number of 
cavities per line. At magnetic fields near H c2 the situation becomes more 
complicated. If the cavity radii are large enough and the distance be- 
tween flux lines is small enough, then it becomes almost impossible to 
assign the pinning sites to particular flux lines. Complicating the problem 
still further is that the flux lines form into bundles. So if a bundle is 
pinned the force on the bundle is the composite Lorentz force on all 
the flux lines in the bundle. Under certain circumstances it may be 
energetically more favorable for the bundle to fragment rather than 
remain pinned. The flexibility of the flux lines also decreases at higher 
magnetic fields because of the repulsive forces between flux lines. 
Therefore, flux lines cannot bend to contact pinning sites. 

The cavity pinning mechanism has also been studied by Silcox and 
Rollins (ref. 86). Silcox and Rollins determined the interaction potential 
between nearest neighbors empirically from magnetization measure- 
ments (refs. 87 and 88) to be 


U(a) 


<Po 


12ttV 3 („ tfc2 „ ) a 2 

\tl c2 — n c \/ 


(173) 


where a is the nearest-neighbor distance between flux lines. Then they 
assumed a density of pinning sites and equated the force of each pinning 
site to Hct<polfar* Using these assumptions they calculated the magnetic 
induction and the associated magnetization curve. Their results showed 
increasing hysteresis with increasing defect content. These results were 
very similar to the experimental results of Livingston and Swartz (refs. 
87 and 88). Narlikar and Dew-Hughes (ref. 89) found that magnetic 
hysteresis in Nb was very sensitive to the dislocation density. Webb 
(ref. 90) calculated the pinning force on an infinitely long flux line due 
to one screw dislocation. He assumed that the nonuniform stress field 
of a dislocation provides the pinning mechanism. He calculated that in 
Nb with a dislocation density of 10 14 per square meter, the pinning force 
would be 2X10 -6 newton per meter. Yamafuji and Irie (ref. 91) have 
criticized Webb’s calculation for assuming static conditions and for 
ignoring the lattice energy of the Abrikosov structure. They analyzed the 
case of a linear flux lattice moving with a constant velocity in the pinning 
potential shown in figure 29. They also assume that the individual flux 
lines in the lattice are displaced from their equilibrium position as the 
line moves through the potential well. The displacement from equilibrium 
is denoted by Ax. If the deformation results from a Hooke’s law force, 
then the equation of motion is 
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MAx + r\Ax + kAx—— 



(174) 


where M is an effective mass per unit length, r\ the viscosity parameter, 
and k a restoring force constant. Considering the energy lost as a flux 
line moves past a pinning site yields a pinning force 

f p =(-kbx) (175) 



As the flux lattice moves past a particular pinning site, it is distorted. 
When the flux line is pulled free of the pinning site the deformation 
energy is transferred to vibrational energy of the crystal lattice. The 
vibrational energy in turn is dissipated to the viscous medium. Recently 
the present status of the relationship between critical current density 
and the pinning site has been reviewed by Livingston (ref. 92). 


PRESENT SITUATION AND UNSOLVED PROBLEMS 

In recent years considerable progress has been made toward a thor- 
ough understanding of the mixed state of superconductors. The present 
understanding is, however, far from complete. The brilliant experiments 
of Trauble and Essmann have proved the existence of the Abrikosov 
lattice. The flux creep theory of Anderson and Kim is in satisfactory 
agreement with the experimental results. Kim, Hempstead, and Strnad 
have correlated the voltage — current — magnetic field relations of type 
II superconductors by experiment and analysis. Their flux motion hy- 
pothesis and the law of corresponding states are in agreement with a 
considerable amount of experimental results. The motion of flux lines 
has been substantiated by the experiments of Hans Meissner. The 
existence of an associated voltage has been affiirmed by Pearl’s experi- 
ments, and Goodman has derived the flux flow voltage from the Ginzburg- 
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Landau equations. The voltage is found to be the voltage that results 
from a Lorentz transformation between a moving system and a stationary 
system. The theoretical analyses of Bardeen-Stephen and Nozieres- 
Vinen have suggested possible mechanisms for the phenomena of flux 
flow. Several mechanisms have been suggested to explain the viscous 
drag postulated by Kim, Hempstead, and Strnad. The results of Tink- 
ham’s oscillating- wave-function damping scheme are extremely close 
to experimental results. 

As impressive as the rate of progress toward a complete understanding 
of the mixed state has been, several areas need further theoretical and 
experimental analysis. The flux flow state must be incorporated into 
either a time-dependent Ginzburg-Landau theory or a still more general 
theory. Attempts toward a time-dependent theory have been made by 
Stephen and Suhl (ref. 93) and by Schmid (ref. 94). Schmid’s solution 
incorporates a relaxation time for the order parameter that is similar to 
Tinkham’s (ref. 83). At the Stanford superconductivity conference John 
Bardeen reported on his attempts to provide general solutions for super- 
conductors where the energy gap parameter varies in space and time 
(ref. 95). 

The experiments of Trauble and Essmann (ref. 48) proved that the 
mixed state is composed of individual flux quanta. More recent experi- 
ments (refs. 96 to 99) utilizing these same techniques have produced 
results which indicate that the triangular lattice of flux lines is not 
always found. Much of this evidence supports an intermediate state in 
type II superconductors. Essmann (ref. 100) interprets these results 
as an indication of an attractive interaction between flux lines. 

Trauble and Essmann (ref. 96) have also studied flux flow in lead foils. 
They applied magnetic field and current and deposited small magnetic 
particles onto the sample. With no transport current through the film 
they found a flux structure consisting of flux bundles of about 50 quanta. 
When a field of 6.0 X 10~ 3 tesla and a current density of 1.2 X 10 8 amperes 
per square meter was applied, the sample was in the flux flow state. 
Under these conditions an evaporation time of about 1 second produced 
the picture shown in figure 30. The dark areas correspond to the regions 
of the sample that the flux lines moved through. The result shown in 
figure 30 is unexpected. If the flux lattice moved through the sample 
as a rigid entity the sample surface would be expected to be blackened 
evenly. Instead, figure 30 shows that the flux lines flow along specific 
paths. Clearly this experimental observation is contrary to the Kim, 
Hempstead, and Strnad flux flow theory, which assumes that all the flux 
lines move together. 

There is other evidence which points to limitations of the flux flow 
model. Hans Meissner has observed that when the current reaches a 
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Figure 30. —Electron microscope picture of flux flow in lead foils (ref. 96). Magnetic field 
intensity, — 6.0 X 10~ 3 tesla; current density, J — 1.2 X 10 8 amperes per square meter. 


certain value he can no longer observe flux moving across the sample. 
This current level corresponds to flux flow velocities equal to the velocity 
of sound. Since superconductivity in the Bardeen-Cooper-Schrieffer 
formulation is a phonon coupling of electrons, it would not be surprising 
at all if the flux lines had some maximum velocity. These ideas were 
presented by Meissner (ref. 57). It has also been observed by Hudson 
and Jirberg (ref. 101) that in niobium the voltage is never directly pro- 
portional to the current (fig. 31). This is contrary to the prediction of 
Kim, Hempstead, and Strnad (eq. (105)). If the number of flux lines in 
motion varies, it is possible to explain almost any nonlinear flux flow 
voltage. Trauble and Essmann have in fact observed instances where 
part of the flux moves while another part is held stationary by pinning. 


Magnetic 

field. 



Figure 31. — Resistive transition near H c2 (ref. 101). 



APPENDIX A 

PLOTS OF CRITICAL FIELD AS FUNCTION OF TEMPERATURE 

The curves below are plots of critical magnetic field as a function of 
temperature for several elements. The source of the curves is reference 
102 . 


HOOxlO ' 4 





APPENDIX 8 

PHOTOGRAPHS OF INTERMEDIATE STATE 

The intermediate state has been optically observed using an experi- 
mental technique that Jakes advantage of the Faraday rotation effect. 
Some pictures taken by Faber (ref. 16) are presented in this appendix. 
In these pictures of an aluminum plate the dark areas are the supercon- 
ducting regions in the materials. 



(a) H — 0 (trapped flux). (b) H — 0.27 H c . 

FIGURE 33.— Intermediate state. 
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(e) H = 0.79 H c . 


(f) = 0.9 He 


Figure 33. — Concluded. 


APPENDIX C 

INTEGRAL TRANSFORMATION 


The purpose of this appendix is to demonstrate the integral transfor- 
mation mentioned on page 32. The integral under discussion is 


IZ=( dr^~ t^ + XKVXfe) 2 ] (Cl) 

J(r>£ ) orr 

Integration can be performed if the integral is first transformed to a 
surface integral. The necessary transformation can be derived from the 
divergence theorem and the vector identity: 

V . (AxB) =B • (VXA) - A * VxB (C2) 

(da- (AxB)=( dr(V ■ (AxB))=( dr(B ■ VxA)- f dr(A - VxB) 


(C3) 


If A = V x h and B = h, then 

£da- ( (V Xh) X h) — J dr(h-VxV + h)-J dr(V X h) - (V X h) 

(C4) 

Recognizing the term with the (V X h) 2 integrand as the second term in 
equation (64) gives 


1Z = 


8 77 


Jdr h 2 — X 2 da- [(Vxh)Xh]+\ 2 J dr[h • (VxVxh)] 


(C5) 


lz — ( J dr{h 2 +\ 2 [h- (VxVx/i)]} — X 2 J da- [(Vxli)Xfi]^ (C6) 


Then recognizing the integrand of the volume integrand as the London 
equation (eq. (6)) yields 


IZ = 


-f 

8 IT J 


du ' [JiX(VxFi)] 


(C7) 
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APPENDIX D 

DERIVATION OF EQUATION (76) 

This appendix shows the steps necessary to derive equation (76) from 
equation (75). In order to get started several equations must be borrowed 
from the Bardeen-Cooper-Schrieffer theory (refs. 26 and 27): 


-H c = 27r i mi 2 A } 


A 



(Dl) 


where N is the density of states at the Fermi energy per unit (energy 
cm 3 ); Substituting for A gives 


H c = 2tt 1 W/ 2 ^ 

Dividing equation (Dl) into equation (75) gives 

i Po_ 

Her 4tt 77 £ X 

H c 27 T^N^fiVF X 2 £ 


Substituting <p 0 = chl2e into equation (D3) yields 


Hci 7 r 1/2 c £ , X 

H c 8 N&ev, X 2 £ 


(D2) 


(D3) 


(D4) 


Then substituting X= (mc 2 l4iTne 2 ) 12 for one of the X’s gives 

Hci _ 7 t n 1 / 2 £ , X 

H c 4 m^Vf Wl* X £ 


(D5) 


FinaDy, replacing N with m 2 vrl27r 2 h 3 , k = X/£ and tv with hlm{iir 2 n ) 13 
can give the desired result: 


He 


77 


H c 


V24 


In 


(D6) 
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Table L— Superconducting Parameters 


Element 

Penetration depth. 

Energy gap, 

Coherence length. 


MO), A 

2A(0) 

io, A 

AI 

500 

0.34 X10- 3 

16000 

Cd 

1300 


7 600 

Hg 

380 to 450 

1.65 


In 

640 

1.05 

4400 

Pb 

390 

2.67 

830 

Sn 

510 

1.15 

2 300 

T1 

920 



V 




Nb 

440 

3.05 

380 

Ta 


1.40 


NbTi 

NbZr 

Nb 3 Sn 


These materials vary considerably 
with composition (refs. 103 and 104) 

i i 
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